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I , Context. Narrow-band surveys to detect Lya emitters are powerful tools for identifying high, and very high, redshift galaxies. 

Q ■ Although samples are increasing at redshifts 2: = 3 — 6, the nature of these galaxies is still poorly known. The number of galaxies 

^ ' detected at redshifts below 2 ~ 3 are also small. 

C/j I Aims. We study the properties of 2 = 2.25 Lya emitters and compare them with those of 2 > 3 Lya emitters. 

. Methods. We present narrow-band imaging made with the MPG/ESO 2.2m telescope and the WFI (Wide Field Imager) detector. 

Using this data, we have searched for emission-line objects. We find 170 candidate typical Lya emitters and 17 candidates that we 
regard as high UV-transmission Lya emitters. We have derived the magnitudes of these objects in 8 photometric bands from u* to 
Ks, and studied whether they have X-ray and/or radio counterparts. 

Results. We demonstrate that there has been significant evolution in the properties of Lya emitters between redshift 2 ~ 3 and z = 
2.25. The spread in spectral energy distributions (SEDs) at the lower redshift is larger and we detect a significant AGN contribution in 
the sample. The distribution of the equivalent widths is narrower than at 2 ~ 3, with only a few candidates with rest-frame equivalent 
width above the predicted limit of 240 A. The star formation rates derived from the Lya emission compared to those derived from the 
UV emission are lower by on average a factor of ~ 1.8, indicative of a significant absorption by dust. 
fSJ ' Conclusions. Lya emitters at redshift z — 2.25 may be more evolved than Lya emitters at higher redshift. The red SEDs imply 

more massive, older and/or dustier galaxies at lower redshift than observed at higher redshifts. The decrease in equivalent widths and 
QQ ' star formation rates indicate more quiescent galaxies, with in general less star formation than in higher redshift galaxies. At 2 = 2.25, 

, AGN appear to be more abundant and also to contribute more to the Lya emitting population. 
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1. Introduction littleprogresshasbeenmade. in practice, the lower redshift limit 

is in fact higher than z w 1.6 because of the drop-off in CCD 

Over the past decade, a large number of so-called Lyman-a sensitivity and a more typical lower redshift limit is z - 2.0. 

(Lya) emitters have been discovered at high redshift. Several ^^^^^ narrow-band surveys have so far been published below 

techniques have been employed in the search for these galax- ^ _ 3. py^^o et al. (1999), Pentericci et al. (2000), Stiavelli et al. 

les, but the by far most common method is that of narrow-band (2OOI), Fynbo et al. (2002, 2003a, 2003b), Francis et al. (2004) 

imaging, where a narrow-band filter is tuned to Lya within a Venemans et al. (2007). Furthermore, Van Breukelen et al. 

particular narrow redshift range. Objects with large equivalent (2005) published a sample of Lya emitters (LAEs) at z ^ 2.5 

widths (EWs) are thus selected by comparing the colours in ^^^^^ integral-field spectroscopy. The difficulty in observing the 

the naii-ow-band image and complementary broad-band images, j^ya line between redshifts 2 < z < 3 lies in the low throughput 

Spectroscopically confirmed Lya emitters now include several of optical systems, the low efficiency of CCDs in this wavelength 

hundreds of sources at redshifts z ^ 3 (e.g., M0ller & WaiTen ^^^^^^ higher extinction in the atmosphere. Even so, the ad- 

1993; Cowie & Hu 1998; Steidel et al. 2000; Fynbo et al. 2001, ^^^^^^^ of a smaller luminosity distance is rewarding in that a 

2003; Matsuda et al. 2005; Venemans et al. 2007; Nilsson et higher ^mx limit equals the same luminosity limit as surveys at 

al. 2007; Ouchi et al. 2008), z - 4.5 (Finkelstein et al. 2007), higher redshift. It also facilitates follow-up observations into the 

z ^ 5.7 (Malhotra et al. 2005; Shimasaku et al. 2006; Tapken nature of these objects 
et al. 2006) and z ~ 6.5 (Taniguchi et al. 2005; Kashikawa et 

al. 2006). However, in the low redshift range, between z « 1.6 At higher redshifts, Lya emitters have been observed to be 

corresponding to the atmospheric cut-off in the UV and z ~ 3, increasingly bluer, younger and smaller with increasing redshift. 

At z ~ 3, Gawiser et al. (2006) inferred stellar masses of a 

Send offprint requests to: knilsson@mpia-hd.mpg.de few X 10® M©, almost no dust extinction, and ages of the order 
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of 100 Myr. In a follow-up paper, Gawiser et al. (2007) stud- 
ied a stacked sample of 52 Lya emitters without Spitzer detec- 
tions and confirmed their results of young ages and that no dust 
appears to be present in these systems, although they inferred 
slightly higher masses. The galaxies in their sample with Spitzer 
detections are presented in Lai et al. (2008). For these galax- 
ies, older ages and higher masses are reported and the authors 
propose that z ^ 3 Lya emitters may have a wide range of prop- 
erties. In Nilsson et al. (2007), a stacked sample of Lya emitters 
at z = 3.15 were studied. Here, small masses and low dust con- 
tents were inferred, but the ages were unconstrained. At even 
higher redshift, Pirzkal et al. (2007) showed that a sample of 
4 < z < 5.7 Lya emitters had very young ages of a few Myr and 
small stellar masses, in the range 10^ — 10® Mq. These results 
agreed well with those of Finkelstein et al. (2007), who studied 
almost 100 Lya emitters at z = 4.5. Finkelstein et al. (2008) 
reported, studying a different sample, finding very dusty, mas- 
sive galaxies with ^1200 as high as 4.5 magnitudes and masses 
as high as several x 10^° Mq, at z — 4.5. They argued that they 
observed a bimodality in the properties of Lya emitters; young 
and blue galaxies versus old and dusty. Thus, the properties of 
z > 3 Lya emitters are poorly constrained, but these galaxies 
tend to be considered young, blue, small and dust-free. With the 
data presented here, we aim to extend the study of the proper- 
ties of Lya emitters to lower redshifts, where a wider range of 
the SED can be studied in the optical/infrared and the luminosity 
distance is smaller, allowing a more detailed analysis. 

This paper is organised as follows. Section |2] presents the 
observations leading to our sample, as well as the data reduc- 
tion. Section[3]includes the object detection and candidate selec- 
tion and a discussion about possible interlopers in the sample. 
In Sect, m we present all the basic characteristics of the sample, 
including photometry, AGN contribution, surface density, and 
sizes, and the equivalent width distribution of the candidates. We 
summarise the results in Sect. |5j 

Throughout this paper, we assume a cosmology with Hq = 
72 km s^i Mpc^i (Freedman et al. 2001), = 0.3 and = 
0.7. Magnitudes are given in the AB system. 



Table 1. Log of imaging observations with WFI. 



Date 


Total exp. 


Average seeing 


10-11.03.2007 


0.56 hours 


1'.'40 


11-12.03.2007 


0.83 hours 


1'.'24 


12-13.03.2007 


1.33 hours 


1'.'04 


13-14.03.2007 


2.50 hours 


1'.'65 


14-15.03.2007 


1.67 hours 


1. 26 


15-16.03.2007 


1.45 hours 


1'.'81 


16-17.03.2007 


3.33 hours 


I'.'ll 


17-18.03.2007 


1.67 hours 


1'.'32 


18-19.03.2007 


3.33 hours 


1'.'15 


19-20.03.2007 


2.50 hours 


1'.'25 


20-21.03.2007 


2.50 hours 


1'.'36 


22-23.03.2007 


1.67 hours 


1'.'13 
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Fig. 1. Transmission of selection filters. The WFI narrow-band 
filter is drawn with a solid line. The dashed line shows the 
SUBARU Bj band filter curve and the dotted fine the CFHT u* 
band filter curve. 



2. Observations 

In March 2007, a 35x34 arcmin^ section of the COSMOS field, 
centred on R.A. = 10''00'"27'' and Dec = 02°12'22'.'7 (J2000), 
was observed with the Wide-Field Imager (WFI; Baade et al. 
1999) on the MPG/ESO 2.2m telescope on La Silla. A log of 
the observations can be found in Table [T] The total dithered 
image consisted of 29 exposures with a total exposure time of 
99624 seconds, or 27.7 hours. The observations were made with 
narrow-band filter N396/12 with a central wavelength of 396.3 
nm and a FWHM of 12.9 nm. This wavelength range corre- 
sponds to z = 2.206 - 2.312 for Lya, z = 0.046 - 0.081 for 
[Oil], and z = 1.52 — 1.59 for CIV and the surveyed comoving 
volumes (after masking, see Sect. 13.11 ) are ^ 329 300 Mpc^"^, 
820 Mpc"^ and 225 000 Mpc"^, respectively. The naiTow- 
band filter curve is shown with the filter curves of the filters used 
for the selection of candidates in Fig.[T] 

The data were reduced using the MPIAphot pipeline devel- 
oped in MIDAS with specific routines to handle WFI data. We 
briefly describe the reduction steps performed on the data. The 
data from the individual CCDs were converted into MIDAS for- 
mat and bias-corrected on each separate CCD. The individual 
CCD frames were then placed into placeholders in an empty 



mosaic, thus creating a full mosaic for each frame including the 
gaps between the CCDs. To correct for bad pixels and columns 
as well as correcting columns with a constant offset compared to 
the surrounding pixels, lamp images with exposure times rang- 
ing between 1 — 220 s were downloaded and analysed. The sci- 
ence frames were then corrected for both the bad columns and 
the offset columns. The frames were then flat-field corrected 
using a master flat created from sky flats taken at the time of 
observations. To determine the offsets between the images, an 
SDSS catalogue of known sources was used, combined with an 
algorithm that identifies sources and matches them with cata- 
logue entries. The information about the shift is then entered 
into the header. In the final step, the images were shifted and 
rebinned onto a gnomonic projection (i.e. by de-projecting all 
great circles onto straight lines), cosmic -ray hits were removed 
and a final, co-added image created. Significant residuals from 
the background subtraction were apparent in the mosaiced im- 
age, in particular around the edges of the individual images used 
to create the mosaic. We used the following procedure to remove 
it. First, we used the "clean" option of the IRAF task imsurfit to 
create a version of the image in which objects had been removed 
and interpolated over We then edited this image interactively to 
remove any residual flux from objects. Subsequently, we fitted 
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a 20x20 piece bicubic spline to this image. We then subtracted 
this fit from the original image. This procedure reduced the back- 
ground variations to a small fraction of the background noise. 
We then flux calibrated the narrow-band image by calculating 
the fluxes of all narrow-band selected objects (see Sect. 13.1b in 
the CFHT u* and the SUBARU Bj images (see Table O and 
interpolated the fluxes in the narrow-band image. We can then 
calculate the zero-point of the image, assuming that the median 
equivalent width of all objects is zero. The 5a detection limit 
in a 3" diameter aperture in the image is 25.3 AB magnitude 
and the 90% completeness limit is 25.1 AB magnitude, corre- 
sponding to Lya luminosities of logL — 42.36 erg s^^ and 
logL = 42.44 erg s^^, respectively, at z = 2.25. 

3. Selection of candidates 

3. 1. Object detection and candidate selection 

For object detection, we used the SExtractor software (Bertin 
& Arnouts 1996). The narrow-band image was used as a de- 
tection image, and objects with a minimum of 8 adjoining pix- 
els and a threshold of 2<t per pixel were selected. For the se- 
lection, we used two broad-band filters corresponding to wave- 
lengths to both the blue and red side of the narrow-band filter, 
the CFHT u* band image and the SUBARU Bj band image, see 
also Fig. [T] Both images were taken from the public data in the 
COSMOS field (Capak et al. 2007). The u* and Bj band im- 
ages were created from combining several smaller sub-images 
to match perfectly the field of the narrow-band image, and they 
were subsequently re-binned and smoothed to match the pixel 
size and source PSF of the narrow-band image. The flux of each 
object detected in the narrow-band image was then measured in 
all three bands in circular apertures with a diameter of 3", where 
each aperture was centred on the position of the object detected 
in the narrow-band image. For the catalogue, we included only 
objects with S/N > 5 in the narrow-band image that were found 
at least 70 pixels from the edge of the image. We also masked 
areas where stray light and bright stars affected the image. This 
left us with a catalogue of 21 275 objects in an effective area of 
^1014 arcmin'^. 

The candidate Lya emitters were selected according to a 
similar method as described in Nilsson et al. (2007). The method 
is based on determining the equivalent width (EW) of a potential 
emission/absorption line located in the narrow-band wavelength 
range. Thus, the flux density was determined in the CFHT u* and 
SUBARU Bj bands and interpolated between the central wave- 
lengths of these filters and the central wavelength of the narrow- 
band filter The EW was then determined by dividing the mea- 
sured flux in the narrow-band with the interpolated flux density. 
Simultaneously, we also calculated the propagated statistical er- 
ror in this measurement. The calculated EW represents a lower 
limit to the true EW of the object as there may be several ob- 
jects, unrelated to the narrow-band object, detected inside the 
aperture in the broad-band images. The EW calculation is also 
complicated by the uncertainties in the continuum interpolation 
(Hayes & Osthn 2006; Hayes et al 2008). To compile the final 
catalogue of Lya emitting candidates, we used two selection cri- 
teria. First, all objects with a measured EW larger than or equal 
to 65 A (corresponding to 20 A in the restframe for Lya) are se- 
lected. To further exclude false candidates at the faint magnitude 
end of the sample, in which the EW measurement is dominated 
by the broad-band flux uncertainties, a significance of the EW 
greater than 2.3 is required. These criteria are designed to detect 
all emission-line objects with large equivalent widths accurately. 
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Fig. 2. Plot of the selection criteria, with narrow-band magni- 
tude on the X-axis and observed equivalent width on the y-axis. 
Selected candidates are shown with error bars, green points cor- 
respond to the GALEX-detected objects, see sec. 13.21 Solid line 
marks the 65 A limit. Objects above this line with EW signifi- 
cance greater than 2.3(7 are selected in the first catalogue. Dots 
mark the entire sample of 5a detections in the narrow-band, and 
dots above the selection criteria without error bars were rejected 
in the visual inspection. 



This final catalogue includes 386 objects. Following the proce- 
dure of Fynbo et al. (2002), three co-authors then individually 
ranked these 386 candidates by visual inspection into three cate- 
gories; first a category of rejected detections (caused by CCD or 
stellar artifacts or enhancement due to airplane/satellite tracks), 
then into two categories of likely candidates and candidates that 
are unlikely to be real but for which no obvious reason was found 
to reject them (see also Fig. |2]l. The 386 candidates were split 
into the different categories of [Likely, Unlikely, Rejected] = 
[187, 61, 138]. In Table |2l we list the coordinates, EWs, sig- 
nificance of the EWs, continuum-subtracted narrow-band mag- 
nitudes Mag396o and mean UV spectral index (see Sect. 14. 11 ) of 
all likely candidates (coordinates of unlikely candidates may be 
obtained from the authors). In the following sections, we only in- 
clude the candidates in the likely category, which provides a total 
of 187 candidates. The Lya magnitudes of our LAE candidates 
are given by the range Magaggo = 22.57 — 25.88 corresponding 
to Lya luminosities in the range of logiLya = 42.13 — 43.45 
(see also the histogram of magnitudes in Fig.[3]l. 

3.2. Interloper discussion 

The only possible contaminants in the sample are [Oil] emit- 
ters at z = 0.06 and CIV emitters at z = 1.56. For [Oil] 
emitters at this redshift, the detection limit corresponds to a 
luminosity limit of logL = 38.7 erg s~^ and the survey 
volume is 820 Mpc~^. For CIV emitters, the corresponding 
limit is logL = 42.0 erg s^^ and the survey volume is 
225 000 Mpc^'^. Unfortunately, no published luminosity func- 
tion for local Universe [Oil] emitters or CIV emitters reach these 
faint luminosity limits, but for [Oil] emitters, extrapolating from 
the luminosity functions of Hogg et al. (1998), Gallego et al. 
(2002), and Teplitz et al. (2003), we conclude that we could 
expect between a few to 200 [Oil] emitters in the survey vol- 
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Table 2. Coordinates, observed equivalent widths, continuum subtracted narrow-band magnitudes and UV spectral index /3 of the 
candidates. Full table is available in the online version. 



Rank 1 candidates 



LAE_COSMOS_# 


R.A. 


Dec. 


EW(A) 


Mag3960 




1 


150.05310 


1.94797 


142.7 ±38.6 


25.50 ± 0.34 


-0.48 ±0.24 


2 


150.21068 


1.95100 


225.8 ±57.0 


25.41 ± 0.32 


-0.82 ±0.33 


3 


150.29437 


1.95107 


564.0 ±43.2 


24.24 ±0.09 


-1.40 ±0.40 


4 


150.25458 


1.95224 


118.4 ±33.4 


25.55 ±0.36 


-2.12 ±0.33 


5 


150.31830 


1.95244 


279.3 ±41.8 


24.91 ±0.18 


-2.08 ±0.44 


6 


150.37505 


1.95619 


95.9 ±24.2 


25.46 ±0.32 


-1.84 ±0.29 


7 


150.15936 


1.95717 


790.3 ± 140.1 


25.04 ±0.21 


-2.97 ± 1.06 


8 


150.37990 


1.96143 


119.6 ±16.8 


24.88 ±0.16 


1.81 ±0.22 G 


9 


150.28378 


1.96167 


248.1 ±50.2 


25.20 ±0.25 


-1.84 ±0.34 


10 


150.30785 


1.96243 


81.0 ±28.5 


25.80 ±0.47 


0.54 ±0.26 



Coordinates are in J2000. Candidates marked with the capital letter G are detected with GALEX, see text below. 
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Fig. 3. Histogram of Lya magnitudes (continuum-subtracted 
naiTow-band magnitudes), Magaggo, of emission-line candidate 
sample. Binsize is 0.2. 



ume. To constrain this value further, the [Oil] luminosity func- 
tion was extrapolated from the Ha luminosity function of Ly et 
al. (2007). Assuming a conversion between [Oil] and Ha lumi- 
nosities given by the star formation rate equations of Kennicutt 
(1998) and Kewley et al. (2004), the [Oil] luminosity Hmit cor- 
responds to an Ha luminosity limit of \ogL — 38.8 erg s^^ in 
this survey. To this limit, Ly et al. (2007) detected 0.2 emitters 
per Mpc^^, corresponding to approximately 160 objects in our 
survey volume, of which only 2% are expected to have an EW 
larger than the selection criteria limit of EWobs = 65 A (Hogg et 
al. 1998). Based on the results of Ly et al. (2007), approximately 
three [Oil] emitters are hence included in the list of candidates. 

We also extracted the photometric redshifts for our candi- 
dates from the catalogue of Gabasch et al. (2008). The catalogue 
is incomplete in areal coverage and covers roughly 80% of our 
field. The catalogue was searched for counterpart objects within 
a 2" radii from the position of our source. If several objects were 
found, the object nearest to our source was chosen as the coun- 
terpart. For a total of 187 candidates, 132 had detected counter- 
parts. Unfortunately, the near-IR coverage of the COSMOS field 
is patchy and shallow, which complicates any photometric red- 



shift determination in the redshift range z ~ 1.3 — 2.5, since in 
this range the Lyman-break has not yet entered the optical win- 
dow and the Balmer break is in the near-IR. This means that we 
are in principle able to exclude [Oil] emitters from the sample, 
whereas other intermediate redshift emitters, such as CIV, are 
difficult to exclude. The redshifts of our 132 detected candidates 
also show that no emitters are consistent with being [Oil] emit- 
ters, and the median redshift of the candidates is z ~ 1.7, with 
most of the candidates being consistent with having a redshift 
of z = 2.25 to within 3ct margins of eiTor To study how many 
CIV emitters may be expected in the survey volume, we used the 
AGN luminosity function of Bongiorno et al. (2007). Since CIV 
emission of restframe equivalent widths larger than EWq > 25 
are only expected to originate in AGN, and with the luminosity 
limits in the optical of this survey, the expected number of CIV 
emitters in our survey volume is '--^ 2. Any contamination of CIV 
emitters in the total sample of candidates is thus minimal. 

The COSMOS field also has publicly available GALEX 
near- and far-UV data. This data may be used to exclude any 
remaining [Oil] emitters from the sample, since the limiting 
magnitude of the GALEX data is ^ 25.5 in the NUV band 
and, assuming a flat spectrum, all [Oil] emitters in the sam- 
ple should be detected in the GALEX data based on the flux 
limit in the Bj band. Lya emitters may or may not be detected 
in the GALEX data, depending on the level of absorption from 
the Lyman forest. Lya emitters would require a transmission of 
less than 60% to be undetected in the GALEX images. M0ller 
& Jakobsen (1990) inferred that approximately 10% of galaxies 
at redshift z — 2.25 have transmissions in the UV higher than 
60%, in agreement with recent observations at z ^ L5 (Siana 
et al. 2007). The conclusion is that we expect all [Oil] emitters 
to be detected in the GALEX images, whereas the total sample 
of GALEX-detected objects is a mix of high transmission Lya 
emitters and [Oil] emitters. Of the 187 candidates, 17 are de- 
tected in the GALEX images (labelled by a G in Table |2]). The 
largest EW of a GALEX-detected candidate is EWobs = 604 A. 
As shown earlier, the expectation is 3 [Oil] emitters in the 
survey volume, and so most of the GALEX-detected objects 
are most likely to be high transmission Lya emitters. We sepa- 
rate the GALEX-detected sample and the non-GALEX-detected 
sample in the following analysis. Thus, we have two samples; 
one robust sample of 170 Lya emitter candidates (removing 17 
GALEX detections from a total of 187 candidates), which is ex- 
pected to be free of interlopers, and one sample of 17 GALEX- 
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Table 3. COSMOS broad-bands used. 



Band 


WUaCl VtlLUi y 


rwT i^A"! 


FWHM ( 


Depth 


— ^ 

u 


CFHT 


3798 


720 


26.4 


Bj 


SUBARU 


4460 


897 


27.3 


g+ 


SUBARU 


4780 


1265 


27.0 


Vj 


SUBARU 


5484 


946 


26.6 


r+ 


SUBARU 


6295 


1382 


26.8 


i+ 


SUBARU 


7641 


1497 


26.2 


Z+ 


SUBARU 


9037 


856 


25.2 


Ks 


KPNO 


21500 


3200 


21.6 



^ The depths are 5o- AB magnitudes as measured in 3" apertures 
(Capak et al. 2007). 



detected objects, possibly including a small contamination of 
[Oil] emitters but consisting mostly of high-transmission Lya 
emitters. The non-GALEX-detected sample is a conservatively 
selected sample defined so that the probability of having non- 
LAE objects in it is minimal, but at the price that it is less com- 
plete. For this reason, it is unsuitable for determining volume 
and surface densities. In Sect. |4.3| we return to this issue and de- 
fine a sample suited for this purpose. From here on, whenever 
the "sample" is mentioned it is referring to the non-GALEX- 
detected sample. 

4. Basic characteristics of LAEs 

4.1. Photometry 

For the photometry of the candidates, we use the public data-set 
from the COSMOS survey (Capak et al. 2007). The photome- 
try presented here has been made in the bands found in Table |3] 
All images are on the common COSMOS tiling grid and were 
cut and re-binned to match the narrow-band image in coverage 
and pixel size. We ran SExtractor in dual image mode, mea- 
suring the flux in apertures defined in the narrow-band image, 
also including the RMS images in the measurement through the 
WEIGHT JMAGE option. The apertures had a diameter of 3". 
Magnitudes were calculated using the zero-point given in the 
COSMOS data release. These magnitudes will, in the case of 
the candidates, be underestimated, since objects unrelated to the 
LAE candidate may be blended into the aperture (see below for 
a discussion of contamination in the aperture). In the following 
we exclude the results in the g+ band as these magnitudes ap- 
pear to be inconsistent. When plotting the SEDs of the objects, 
the g+ band magnitudes are consistently lower than the other 
SED points by 0.1 — 0.5 magnitudes. 

We calculated the UV spectral index of (3 for the candidates, 
which is defined to be: 

h « (1) 

where f\ is the flux density per wavelength interval in the UV 
range. The spectral index is calculated as in Meurer et al. (1997) 
and Hathi et al. (2008) for the Bj, Vj, and r+ bands using the 
following equations: 

/3Bj-Vj = 4.456 x (Bj-Vj)-2 (2) 
/3Bj-.+ = 2.673 X (Bj - r+) - 2 (3) 

where Bj, Vj, and r+ are the AB magnitudes in the respective 
bands. A plot of the results can be seen in Fig. |4] and the to- 
tal (3, calculated to be the mean of the two measurements as 




Fig. 4. UV spectral index calculated from the Bj— Vj and Bj— r+ 
colours. The two colours give approximately the same results, 
as seen by the solid line marking a linear relation. The spread in 
the spectral index is very large. AGN (see sec. 14.2b are marked 
in red. Objects with red SEDs can be found in the upper right 
corner and objects with blue SEDs in the lower left corner. 



described in Eq. |2] and |3] for our LAE candidates is found in 
Table |2] Most of the candidates are within the expected range of 
— 2 < /3 < 0, with some outliers at both larger and smaller (3. 
Of the four objects with Psj-vj > 2, one is also detected in the 
Ks images, and two are affected by contamination from nearby 
objects. These objects are in the mid range of narrow-band mag- 
nitude, but in general have higher fluxes in the Vj band. On the 
other end of the scale, the three objects with jSnj-Vj < —4 are 
all drawn from the faint end of the magnitude distribution in both 
the narrow-band and the Vj band, which explains their large er- 
ror bars. 

The SEDs of the candidates have characteristics that are 
atypical of higher redshift LAE SEDs. It has previously been 
reported that LAEs have SEDs that are in many cases blue 
(Gawiser et al. 2006, 2007; Pirzkal et al. 2007). Only two 
other surveys reported detections of red SEDs in several LAEs 
(Stiavelli et al. 2001; Lai et al. 2008), although one extremely 
red object was found by Nilsson et al. (2007). Among the 170 
objects detected in this survey, a total of 118 candidates have 
Vj— i+ colours greater than zero, indicative of red colours, and 
five candidates have undefined colours due to upper limits in 
both bands, or only in the i+ band. The median values of Bj— Vj 
and Vj— i+ among candidates with 3(t detections in both bands 
is 0.14 and 0.15, respectively. In Fig.|5] the colours Bj— Vj and 
Vj— i+ are plotted in a colour-colour diagram. In this figure, the 
tracks of a single stellar population, with increasing age and dust 
content, are superimposed. These populations were simulated 
with GALAXEV (Bruzual & Chariot 2003) in age steps of 1, 10, 
50, 100, 500 and 750 Myrs and 1 and 2 Gyrs. The dust extinction 
was increased in 9 steps of E{B — V) between 0.0 and 1.0 with 
increments of 0.125 using the Calzetti extinction law (Calzetti 
et al. 2000). We also completed simulations of an exponentially 
decaying star formation rate (SFR) and constant SFR. However, 
the colours of these populations occupy a smaller region of the 
diagram than those of the single stellar populations and were 
omitted from the plot. Most of the candidates are consistent with 
being young galaxies with ages in the range 10^ — 10^ yrs with 
various dust contents. Some objects appear evolved with a range 



6 



K.K. Nilsson et al.: Lya emitters at z = 2.25 



2.0 

1.5 - 
1.0 - 



I 0.5 : 

> : 

0.0 7 
-0.5 7 

-1.0 L_ 

-1.0 



I , ^ 5.0, 




10" 



10" 



10" 




log(X„„) (/2m 



Fig. 5. Colour-colour diagram with Bj— Vj on the x-axis and 
Vj— i+ on the y-axis. Objects with upper limits (on the flux) in 
one band are shown as arrows with 3a upper limits and objects 
with a colour in either bands that is 1.5a away from the me- 
dian colour in that band are shown with error bars. AGN are 
marked in red. The solid lines display the evolution of a single 
stellar population with increasing dust and age; age causes the 
objects to move predominantly to the right in the plot and dust 
predominantly upwards. The markers 5.0 and 7.5 correspond to 
age steps of 5.0 x 10* and 7.5 x 10* yrs and the dust steps go 
from E{B - V) of 0.0 (point a in the plot) to 1.0 (point b in the 
plot) in steps of 0.125. 



Fig. 6. Stacked SED of the sample. Points are the stacked mag- 
nitudes in six bands {u*, Bj, Vj, i+, z+ and Kg) and in three 
bins; black points are the full sample, excluding i^s-detected 
galaxies and galaxies with bright nearby objects, red stars are 
those with red Vj— i+ colours and blue points those with blue 



Vj- 



colours. Red diamonds mark the stacked i^, -detected 



sample without AGN-detected galaxies. All are normalised to 
the Bj flux of the Nilsson et al. (2007) SED. The lines are two 
best-fit SEDs of z ~ 3 LAEs from Nilsson et al. (2007; green 
fine) and Gawiser et al. (2007; black fine). 



Table 4. Stacked magnitudes for the sample. 



of dust content. There are also some objects that appear evolved 
but fall below even the zero-dust line in the plot, which may be 
affected by nearby objects. Another group of objects are very 
red in terms of their Vj— i+ colours but are too blue in Bj— Vj 
to be explained by simple stellar populations. These may be re- 
produced more accurately by models with more complicated star 
formation histories. 

In the sample, 12 objects have A'^ band detections above the 
limiting magnitude of 22.2 {3a, AB), corresponding to 7% of 
the complete sample. Note that four of these 12 A's -detected ob- 
jects are also selected as AGN (see Sect. 14.21 ). The ATg-detected 
objects are all drawn from the brighter end of the candidate sam- 
ple, but have a range in Vj magnitude of 21.5 — 26, whereas the 
entire candidate sample have a range in magnitude of 21.5 — 27, 
i.e. not only the brightest candidates have Kg detections. We 
stacked the thumb-nail images of the candidates without Kg 
detections, removing objects from the stack with nearby, unre- 
lated, bright detections. The total stack consisted of 144 can- 
didates (with GALEX detections and objects with nearby, un- 
related detections removed) and revealed a mean detection in 
the sample of magnitude 24.00 ± 0.29. The sample was also 
divided into two subsamples with Vj— colours Vj— i+ > 
and Vj— 1+ < 0, respectively. In the subsample with red Vj — i+ 
colours, consisting of 96 of the 144 candidates, a clear detec- 
tion with magnitude 23.65 ± 0.19 was made. In the subsample 
with bluer colours, no detection was made and an upper limit of 
24.26 (3a) was derived. In Fig.|6l we show the stacked magni- 
tudes of the sample divided into four bins; the total sample, the 
sample with Vj— > (called the "red" sample), the sample 
with Vj— i+ < (called the "blue" sample), and the ATg-detected 
sample. The total sample consists of 144 candidates, the red sub- 





Total sample 


Red sample 


Blue sample 


Jfs -detected 


u* 


25.21 ± 0.03 


25.16 ±0.05 


25.34 ±0.06 


24.75 ± 0.06 


Bj 


25.14 ±0.03 


25.11 ±0.04 


25.34 ±0.05 


24.00 ± 0.02 


Vj 


24.99 ± 0.08 


24.89 ±0.11 


25.18 ±0.03 


23.04 ±0.01 


i+ 


24.95 ±0.14 


24.72 ±0.13 


25.37 ±0.15 


21.83 ±0.01 


z+ 


24.86 ±0.22 


24.57 ±0.15 


25.40 ±0.16 


21.46 ±0.01 


Ks 


24.00 ± 0.29 


23.65 ±0.19 


25.55 ± 0.98* 


20.44 ± 0.02 



Upper limits are 3(7. -detected sample is without AGN-detected 
galaxies. *BIue sample Ks magnitude is calculated from the Kg 
magnitudes of the total and red samples. 



sample consists of 96 candidates, and the blue subsample con- 
sists of 48 candidates. The magnitudes in Fig. |6] are also found 
in Table |4] In Fig.|6] it is seen that the stacked SED of the total 
sample is redder than those at z 3; the points of the Tion-Kg- 
detected galaxies are in principle consistent with the Nilsson et 
al. (2007) fit within 2a but inconsistent across the spectrum with 
the Gawiser et al. (2007) fit. The normalised Kg magnitude for 
the total sample is 25.45 ± 0.29, whereas the same magnitudes 
measured from the synthetic SEDs of Nilsson et al. (2007) and 
Gawiser et al. (2007) are 26.00 and 26.39, respectively. Hence, 
roughly two-thirds of the entire sample exhibit on average more 
red SEDs than at z ~ 3. Noteworthy is also that SED of the blue 
sample (with Vj— ■«+ < 0) is consistent with the z ~ 3 SEDs. 
The N]~Kg slopes of the red and blue samples are 1.24 ± 0.22 
and —0.37 ± 0.98, indicating a clear difference in the properties 
of the blue and red samples of galaxies. For reference, the VyKg 
slope of the ATg-detected sample is 2.68 ± 0.04. A large spread 
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in colours is evident in this survey, from the "classical" blue, flat 
SED galaxies to galaxies with very red SEDs, where the redder 
SEDs constitute two-thirds of the sample (see also Finkelstein et 
al. 2008 for a similar result). In a future paper, the Spitzer fluxes 
of these galaxies will be studied. 

To study the level of contamination in the apertures, and the 
possibility that the very red SEDs originate from this source, the 
narrow-band catalogue was searched for other detected objects 
within two aperture radii (i.e. within 3") of each LAE candi- 
date. Of the 170 non-GALEX-detected candidates, 26 have one 
or more detected objects within this search radius. This corre- 
sponds to a potential contamination rate of 15 % in the photo- 
metric measurements. Two objects have more than one nearby 
object. We also wish to determine if the apparent high number 
of red objects is caused by flux contamination from unrelated ob- 
jects in the aperture. Of the 26 candidates with nearby objects, 
19 have a Vj — i+ colour in the aperture greater than zero, and 15 
have Vj— i+ colours larger than the median of Vj— i+ for all 170 
candidates. This implies that the contamination is not confined 
to only apparent red objects but affects all types of SEDs. We are 
thus confident that the contamination in the photometry is small, 
and that contamination cannot explain the apparent red colours 
of the LAE candidates. 

4.2. AGN contribution 

The COSMOS public data includes a release of X-ray data in 
three bands, 0.5 - 2.0, 2.0 - 4.5, and 4.5 - 10.0 keV, as ob- 
served by the XMM observatory (Hasinger et al. 2007). We used 
the catalogue of objects of Cappelluti et al. (2007) to measure 
the AGN fraction of all candidates. The catalogue consists of 
1390 entries representing point-like sources and gives the fluxes 
in the three bands including errors. The catalogue detects objects 
to limiting fluxes of 7.2 x lO^^^, 4.0 x 10"l^ and 9.7 x lO^^^ 
erg s^^ cm^^, respectively, in the bands 0.5 — 2.0, 2.0 — 4.5, 
and 4.5 — 10.0 keV with a confidence of roughly 4.5f7 (Hasinger 
et al. 2007). This catalogue was searched for objects detected 
within a radius of 8" (i.e. 2 pixels in the XMM images) of the 
candidates. Fifteen of these potential counterparts were detected, 
of which five were from the list of GALEX-detected candidates. 
The X-ray fluxes of these candidates are found in Table |5] Note 
that the probability of finding an X-ray object within an aper- 
ture of radius 8" is ~ 0.009, resulting in potentially two random 
detections in the combined GALEX and non-GALEX-detected 
sample. For the objects with detections in two or all bands, the 
hardness ratios can be calculated. Hardness ratios are defined to 
equal: 

HRl = (S2.0-4.5 — So. 5-2.0) / (S2.0-4.5 + So. 5-2.0) (4) 

HR2 = (S4. 5-10.0 ^ S2.0-4.5) / (S4. 5-10.0 + S2.0-4.5) (5) 

The hardness ratios, when applicable, can be found in Table |5] 
All measured hardness ratios (HRl) are positive, indicating that 
the LAE candidates are type 2 AGN (e.g. Norman et al. 2004). 
The low luminosity of the GALEX-detected candidates, assum- 
ing that they are [Oil] emitters indicate that they are normal, 
star forming galaxies, or that they are Lya emitters. The total 
sample of GALEX-detected and non-detected candidates with- 
out individual X-ray detections were stacked, and produced no 
detection in any band in a stack of 170 objects (the nine indi- 
vidual detections were removed from the stack, as well as eight 
candidates with unrelated objects detected in the vicinity). The 
non-detection in a stack of 170 objects imply an upper limit to 



the mean X-ray flux of approximately a factor of ten smaller than 
those in the catalogue of Cappelluti et al. (2007). 

We also searched the catalogue of Schinnerer et al. (2007) 
for radio counterparts to the candidates. The data on which the 
catalogue is based was taken with the VLA array and covers the 
entire COSMOS field to a depth of a few x lO^Jy. The cata- 
logue presented in Schinnerer et al. (2007) contains 3643 source 
detections. This catalogue was searched with the same criteria 
as for the X-ray catalogue (i.e. detections within 8" of the candi- 
date positions) and six objects were found. The number of ran- 
dom detections in a 8" radius aperture is 0.026. Of the six radio 
detections, one object is GALEX-detected (candidate 173). The 
radio fluxes of these sources can be found in Table |5] Only one 
source has a detection in both the X-ray and radio data. This is 
in principle not a problem since AGN can be both radio-loud 
and radio-quiet. The radio objects without X-ray detections are 
more enigmatic, but can be explained by the relatively shallow 
luminosity limit of the X-ray observations. If the radio fluxes 
are converted into star formation rates using the conversion rate 
of Condon (1992), they correspond to values in the range of 
~ 400— ^ 1300 Mq yr^^ whereas the upper limit to star for- 
mation rates in the XMM survey is ^ 1200 Mq yr^^, using the 
conversion of Ranalli, Comastri & Setti (2003). 

The AGN fraction of the non-GALEX-detected LAE sample 
is thus at least ^ 5%, if it is assumed that the X-ray and radio ob- 
jects within 5" of the narrow-band sources are true detections. It 
is clear that the shallow flux limit to the deepest X-ray band, cor- 
responding to a luminosity of Zog(Lo. 5-2.0) = 43.4 erg s~^, en- 
sure many AGN of lower luminosities to be undetected, as con- 
firmed by the radio detections without X-ray counterparts. The 
AGN fraction should therefore almost certainly be even higher, 
although we did not detect any X-ray flux in the mean, stacked 
total GALEX and non-GALEX-detected sample of 170 objects. 
Previous studies of LAEs at higher redshifts exhibited very small 
AGN contributions, from less than one percent in the z ~ 4.5 
LAEs studied by Wang et al. (2004) to 1 % at z 3 as in- 
ferred by Gawiser et al. (2007). Lehmer et al. (2008) calculated 
the 'AGN fraction function" of Lya emitters at z ~ 3 both in the 
field and the overdense region of SSA22 (Fig. 3b in their paper). 
If we convert our 0.5 — 2.0 keV detection limit to their 8 — 32 keV 
luminosity range, we should expect to detect an AGN fraction of 
0.5 — 1.5 %, where the lower number is related to their field 
result and the higher to the SSA22 result. Based on this compre- 
hensive study of X-ray detected Lya emitting AGN, it is clear 
that the AGN fraction detected here is larger than at z 3, as ex- 
pected from the Lehmer et al. (2008) survey. Ouchi et al. (2008) 
argued that the missed fraction in AGN X-ray searches among 
LAEs may be as high as 10%. We performed the same test as in 
the Ouchi et al. (2008) publication. Based on the quasar SEDs of 
Elvis et al. (1994) and the UV spectrum of quasars of Telfer et 
al. (2002) and Richards et al. (2003) we inferred the ratio of ob- 
served 0.5 — 2.0keV flux to Lya flux of 1.87 and 1.51 for radio- 
quiet and radio-loud quasars, respectively . For the flux limit in 
the X-ray band, the maximum magnitude observed for Lya is 
then 23.34. In our complete GALEX and non-GALEX-detected 
samples, nine objects have magnitudes brighter than this, four of 
which are GALEX-detected. Of the nine objects, six have X-ray 
detections (all four GALEX-detected sources and two LAE can- 
didates). Thus, if we exclude the four GALEX-detected sources, 
two out of five Lya bright galaxies host AGN and the fraction 
of X-ray detected AGN in the non-GALEX-detected sample is 
40%, four times the fraction found by Ouchi et al. (2008) for 
z ^ 3 LAEs. To conclude, our detection of a 5% AGN contribu- 
tion in the non-GALEX-detected LAE sample is consistent with 
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Table 5. X-ray and radio properties of counterparts of the candidates. 



LAE candidates 



LAE. 


So. 5-2.0 


S2.0-4.5 


S4. 5-10.0 


HRl 


HR2 


log(Lo. 5-2.0) 


Offset 


Sl.4GHz 


Offset 


COSMOS.* 














(X-ray) 




(Radio) 


25 


129.0 ± 1.87 


191.0 ±5.63 


133.0 ± 7.03 


0.19 ±0.01 


-0.18 ±0.03 


45.67 ±0.02 


11.06 





— 


36 


— 


— 


— 


— 


— 


— 


— 


140 ± 27 


0.10 


57 


— 


— 


— 


— 


— 


— 


— 


47± 10 


6.17 


82 


1.25 ± 1.97 


— 


— 


< 0.52 


— 


43.66 ± 0.41 


1.07 


— 


— 


101 
















137 ± 27 


0.06 


113 


1.21 ±4.10 


9.58 ±29.1 




0.78 ±0.61 


< 0.01 


43.64 ± 0.64 


10.55 






115 


7.20 ±6.00 


16.2 ± 21.5 




0.38 ± 0.59 


< -0.25 


44.42 ± 0.78 


0.17 






119 


0.92 ± 0.26 






< 0.63 




43.55 ± 0.14 


2.67 






121 


17.6 ±0.73 


31.0 ± 2.54 


18.3 ±3.52 


0.28 ± 0.04 


-0.26 ±0.11 


44.81 ± 0.03 


6.36 






129 


1.10 ±0.27 






< 0.57 




43.60 ± 0.12 


1.95 






132 
















47± 10 


9.83 


140 


17.5 ±0.72 


46.3 ± 2.97 


44.9 ± 4.73 


0.45 ± 0.03 


-0.02 ±0.06 


44.80 ± 0.02 


0.74 






150 


2.6 ±0.30 


8.3 ± 1.48 




0.52 ±0.07 


< 0.08 


43.98 ± 0.05 


0.62 


107 ± 30 


0.21 


161 


2.42 ±0.33 






< 0.25 




43.95 ± 0.06 


0.21 






GALEX-detected candidates 


42 


4.9 ±2.92 






< -0.10 




40.6 ± 0.39 


0.55 






108 


8.3 ±0.60 


16.9 ± 2.14 




0.34 ±0.06 


< -0.27 


40.83 ± 0.03 


1.01 






126 


7.9 ±0.50 


12.0 ± 1.79 




0.38 ±0.59 


< -0.11 


40.81 ± 0.03 


0.45 






155 


3.0 ±0.34 


5.6 ± 1.21 




0.30 ±0.10 


< 0.27 


40.39 ± 0.06 


1.57 






173 
















110 ±25 


9.44 



Data from the catalogues of Cappelluti et al. (2007) and Schinnerer et al. (2007). X-ray fluxes are given in units of 10~^^ erg s~^ cm~'^ 
and radio fluxes in units of /ijy. HRl and HR2 correspond to the hardness ratios defined in Eq.|4]and[5] The X-ray luminosity is in erg s~^, 
assuming the objects are located at redshifts z — 2.25 and 0.06 for LAE and GALEX-detected candidates respectively. Offsets are distance 
from centroid of narrow-band detection to X-ray/radio detection in units of arcseconds. It is expected that the objects with offsets larger than 
a few arcseconds (25, 57, 113, 121, 132 and 173) are unrelated to the candidates due to the large offset and the relatively large probability of 
random association. 



previous results, but is indicative of a higher AGN fraction being 
present at this redshift. 

4.3. Surface density, sizes and SFR 

For the surface density determination, we included both the 
GALEX-detected and non-detected objects to obtain as com- 
plete as possible a measure of the number density. Since our 
image has areas of poorer signal-to-noise, as well as stellar arti- 
facts, we selected sub-images of superior quality with a total area 
of ~ 288 arcmin'^, corresponding to 28% of the total surveyed 
area, to use for the surface and volume density calculations. In 
these areas, the selection is complete, and we found 54 candi- 
dates within these sub-images. This implies a surface density of 
LAE candidates of 0.19 arcmin^^ Az^^, or 1.91 arcmin^'^ z^^. 
The volume density is 0.00062 Mpc^^^, which is in the lower 
range of that observed at redshift z ^ 3. Fynbo et al. (2001) 
summarise the surface densities of several early z ^ 3 sur- 
veys. With the exception of the Steidel et al. (2000) survey in 
the overdense SSA22 field, these surveys all determined values 
of 2.11 — 5.9 arcmin"^ z^^ to the flux limit of this survey, but 
with large error bars. Hence, the values at z > 3 are consis- 
tently higher than our value, but in most cases the measurements 
agree to within la. For LAE candidates brighter than the (5a) 
luminosity limit of this survey, Nilsson et al. (2007) identified 
six candidates in their survey, coiTesponding to 0.0018 Mpc^^^, 
which represents a decrease of roughly a factor of three although 
based on a small sample. Gronwall et al. (2007) determined a 
space density of 0.00057 Mpc^^^ above the 5a luminosity limit 
of our survey, which is consistent with the space density found 



here. Stiavelli et al. (2001) presented a survey of z = 2.4 LAEs 
and found a volume density of 0.0001 Mpc^'^ above a luminos- 
ity limit of logi = 42.93 erg s~^. The same number for this 
survey is 0.00009 Mpc^'^. The numbers agree reasonably well 
with previous results. Finally, Prescott et al. (2008) completed 
a survey of z ~ 2.7 LAEs around a so-called Lya blob (e.g. 
Steidel et al. 2000; Matsuda et al. 2004; Dey et al. 2005; Nilsson 
et al. 2006). They argued that the central part of their field is 
overdense, and determined a number density of 0.0021 Mpc^'^. 
This is roughly three and half times higher than in our survey. At 
the edge of their field, the density is instead 0.0012 Mpc^'^, in 
closer agreement with, but still higher than, our result. Thus, the 
surface density in this survey is in almost all cases lower than 
in higher redshift surveys, but could also be consistent with pre- 
vious results. If one considers that number densities have been 
found to vary by factors of 2 — 5 in 0.2 deg^ fields at redshift 
z = 3 (Ouchi et al. 2008), it is clear that more data at both 
higher and lower redshifts is needed to resolve this issue finally. 

As in Nilsson et al. (2007), the sizes of the candidate Lya 
and GALEX-detected objects measured in the naiTow-band and 
r+ broad-band images are presented in Fig. |7] Six LAE ob- 
jects are excluded due to non-detections in the r+ band. The 
FWHM was calculated using the FLUXJIADIUS output from 
SExtractor, and the PSF in the narrow-band image was calcu- 
lated by using 9 objects in the image with fluxes in the range 
of the candidates and the SExtractor parameter CLASS_STAR 
larger than 0.8, which inferred a PSF of 0.96". The PSF of the 
r+ image is similar to that of the narrow-band image PSF (the 
PSF of the r"*" image is 1.06"). As can be seen in the figure, 
we have two marginal Lya blob detections in this survey (of 
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Fig. 7. Size distribution of candidates as measured in the narrow-band (hatched histogram) and r+ (empty histogram) images for 
LAE candidates and GALEX-detected candidates respectively. Sizes as measured by the FLUX_RADIUS command in SExtractor 
The solid line represents the PSF of the narrow-band image, with the dotted lines showing the la of the PSF measurement. These 
values are very similar to those of the r+ image, see text. 



FWHM 3.46" and 3.70"), if a Lya blob is defined to have a ra- 
dius larger than 15 kpc (corresponding to a diameter of 3.75"). 
These two (candidates 54 and 66) are also extended in the broad- 
band images with FWHM ^ 2" and are thus "normal" extended 
Lya emitters. We conclude that we have no blobs in this survey. 
Using the results from Matsuda et al. (2004), Saito et al (2006), 
Gronwall et al. (2007), and Prescott et al. (2008), the expected 
number density of Lyman-alpha blobs ranges from a few to a 
few hundred x 10^^ Mpc^'^. The upper limit to the survey is 
3 X 10^^ Mpc^^, and hence at the low end of the predic- 
tions from previous results. Matsuda et al. (2004) and Prescott et 
al. (2008) argued that Lya blobs reside in overdense regions of 
space. If this assumption is correct, it agrees well with the two re- 
sults from this survey that i) we detect a tentatively lower number 
density of LAEs at this redshift than expected from z ~ 3 obser- 
vations, and //) no Lya blobs are found in this survey. Both these 
observations indicate that the COSMOS field is under-dense at 
z — 2.25. We also note that there is a rapid decline in FWHM 
at larger radii in Fig. [T] indicating that Lya blobs and LAEs are 
two separate categories of Lya emitting objects. If this was not 
the case, a tail of objects should be seen towards larger radii. No 
sample has yet been sufficiently large to confirm whether Lya 
blobs and LAEs are separate populations. It is also apparent that 
the candidate emitters are in general always more extended in the 
narrow-band image than their broad-band counterparts, i.e. the 
broad-band counterparts are consistent with being pure point- 
sources whereas the narrow-band objects are significantly more 
extended than point-sources (cf. M0ller & Warren 1998; Fynbo 
et al. 2001, 2003). This may arise from diffuse scattering of Lya 
within the galaxy (see also Laursen & Sommer-Larsen 2007) 
and/or due to the diffuse ISM (Ostlin et al. 2008). The extent 
of this scattering is expected to be much smaller for [Oil] emis- 
sion, hence, the detection of significantly extended emission in a 
GALEX-detected object may therefore be indicative of it being 
a Lya emitter. 

The SFRs of the candidates are calculated using the com- 
mon star formation indicators of the rest-frame UV continuum 
emission at A = 1500 A and the Lya emission. The Lya emis- 
sion used in this case was the continuum-subtracted narrow -band 



flux. The conversion rates that we use are the following: 



or -ttuv — 



/ . X 47rd? X 1.4 • 10-28 



1 + z 



-Moyr-i (6) 



where the flux density f = is in erg s ^ cm ^ Hz ^, dr is the 
luminosity distance in centimeters and z the redshift (Kennicutt 
1998). For Lya, the conversion is: 



SFR 



'Lya 



8.7x1.12-1041 "-^ 



(7) 



where ^Lya is the flux in the Lya line in erg s~i cm~^ 
(Brocklehurst, 1971, Kennicutt 1983). A plot of the two result- 
ing SFRs can be seen in Fig. [8] The SFRs are in the range of 
^1 — 30 Mq yr~i in the Lya measurement, some of which are 
unusually high for what is generally seen in Lya emitters. The 
dashed line in Fig. [8] corresponds to a 1:1 relation between the 
UV/Lya derived SFRs. A correlation is seen between the two 
measurements. The median of the ratio of UV-derived SFR to 
Lya-derived SFR is 1.8. Looking at Eqs.|6]and|7] it is clear that 
a relation exists between the ratio of the two SFR measurements 
and the equivalent width of an object, such that: 



SFRuv cons. 

oc 



SFR 



■Lya 



EWn 



(8) 



In this equation, E'VFo is the rest-frame equivalent width of an 
object. The constant in the nominator depends on the spectral 
slope of the object, and on dust properties within the galaxy. For 
a flat slope in and no dust this constant equals 67 A. For this 
assumption, our median ratio of 1.8 corresponds to a median 
EW of 37 A, which is consistent with the measured median EW 
of 41.6 A. The ratio is higher than that found by Ouchi et al. 
(2008), who measured a typical ratio of ~ 1.2 in a sample of 
z = 'i.l LAEs, and Gronwall et al. (2007), who, after correct- 
ing for the (1 + z) factor, also obtain a typical ratio of ^ 1.0 
(Gronwall, priv. communication). These results are also related 
to those presented in Sect. 14.41 We note that several objects have 
ratios of below one. In total, 19 objects have measurements that 
deviate by up to ia from the 1:1 line. The explanation for these 
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Fig. 8. Star formation rates as derived from the UV rest-frame 
1500 A emission and the Lya emission for the sample of the 170 
candidates. Red points indicate the AGN in the sample, plotted 
for reference. The measurements correlate well. The solid line 
indicates the median of the fraction of SFR from the UV and Lya 
in the sample, and the dashed line indicates a 1 : 1 ratio between 
the two measurements. 



measurements may originate in the age dependence of the UV 
SFR conversion factor According to Kennicutt et al. (1998), the 
constant in Eq.|6]is valid for galaxies with constant star forma- 
tion in the past > 10^ yrs. For starbursts of younger ages, the 
conversion factor should increase. This would bring the ratios 
closer to the 1 : 1 line. It appears true for at least 1 1 % of the sam- 
ple (19 out of 170), which show signs of young starburst activity. 

4.4. Equivalent width distribution 

At higher redshift (z > 4), emitters with EWs larger than 240 A 
have been predicted (Schaerer 2003; Tumlison et al. 2003) and 
observed (Malhotra & Rhoads 2002; Shimasaku et al. 2006; 
Stanway et al. 2007). These large EWs supposedly cannot be ex- 
plained by a normal initial mass function but have been proposed 
to be signs of Population 111 stars. At redshift z ^ 3, small sub- 
samples of detected emitters exhibit extremely high EWs (a few 
percent in the surveys of both Gronwall et al. 2007 and Nilsson 
et al. 2007, and up to 20% in Ouchi et al. 2008). In Fig.lH the dis- 
tribution of rest-frame equivalent widths of the candidates can be 
seen, also divided into the red and blue sub-samples. The sam- 
ple includes six emitters with restframe EWs larger than 240 A, 
corresponding to ^ 4%. 

In Gronwall et al. (2007), the data is reproduced with an ex- 
ponentially declining function: 

N = Cx e-i^w/^o (9) 

where N is the number of objects, C is a normalisation constant 
and Wo is the e-folding length of the distribution. Gronwall et al. 
(2007) inferred that wq ^ 76t" A. The EW distribution pre- 
sented here was fitted here by a best-fit solution minimising the 
X^, and the errors of the fit were calculated by producing 100 000 
distributions, normally distributed within the margins of errors in 
the equivalent widths. Each cumulative distribution was repro- 
duced by an exponential function. No correction for filter trans- 
mission curve was done since the filter is nearly box-shaped (see 
also Fig.[T]). The best-fit solution to the e-folding EW for the total 



non-GALEX-detected sample is uiQ. tot = 48. 5± 1.7 A, less than 
two thirds of that inferred by Gronwall et al. (2007). This implies 
that, at z — 2.25, there is a higher fraction of LAEs with small 
EWs than at z = 3.1. Corresponding values for the red and blue 
sub-samples are wo.red = 42.2 ± 1.7 and wqmub = 61.5 ± 4.5, 
respectively. The fits indicate that the redder galaxies have a nar- 
rower distribution than the bluer galaxies, a result also consistent 
with the redder galaxies being dustier than the bluer galaxies. 
This relation was previously observed for LBGs (Shapley et al. 
2003). 

Since the Lya EW of a starburst declines with time (Chariot 
& Fall 1993; Schaerer 2003) one might seek to interpret the 
change in width of the distributions in the redshift interval z — 
3 — 2 as a simple evolution of individual objects. To test this sim- 
ple idea, we attempted to recreate both the z = 3.1 and z = 2.25 
EW distribution from the simple evolution of single stellar pop- 
ulation starbursts by using both the predictions of Chariot & Fall 
(1993) and Schaerer (2003). We found that neither of the two 
distributions can be recreated in this way. The decline in the 
Lya EW of a starburst is so rapid that a cumulative EW dis- 
tribution similar to those in Fig. |9] will be much narrower than 
those observed at both redshifts. The only way of obtaining dis- 
tributions similar to those observed is to allow the majority of 
the galaxies to have continuous star formation. We can therefore 
conclude that at maximum a small fraction of the LAEs at red- 
shifts 2 = 2 — 3 are single starburst objects, and are most likely 
better explained by alternative star formation histories. Finally, 
the shift in the width of the distribution between z = 3 and 
z = 2 cannot be explained by simple evolution with time, but 
is further evidence of a higher dust content in LAEs at redshift 
z ~ 2. Worth noting is that a narrower EW distribution is con- 
sistent with the higher UV-to-Lya SFRs described in Sect. 14.31 
since the EW is proportional to the ratio of the flux in the Lya to 
the flux density in the UV, similar to the SFRs. 

5. Conclusion 

In this paper, the observations and selection leading to a sample 
of 170 robust LAE candidates at z = 2.25 have been presented. 
This is the first large-scale survey of LAEs at redshifts around 
two, and several indications of evolution are seen in the proper- 
ties of LAEs between redshift z 3 and z ~ 2. This is however 
unsurprising, since the separation in these redshifts corresponds 
to longer than 1 Gyr The main arguments for an evolution in the 
properties of this type of galaxy are; 

- The large spread in SEDs 

The most conspicuous evolution detected was that of the emer- 
gence of a large number of "red" LAEs. The 12 individual detec- 
tions in the relatively shallow Ks images and the red mean SED 
of the total stacked sample, presented in Fig. |6] indicated that 
a significant subsample, of up to two-thirds of the total sample, 
of LAEs are massive, and potentially dusty, galaxies. We also 
showed that the rest-frame UV colours of these galaxies cannot 
be explained by a simple stellar population, and that even though 
the majority of the candidates are consistent with being young 
galaxies, a significant fraction of the sample displays more com- 
plicated colours. 

- Potentially larger fraction of AGN 

The result above is also related to the detection of at least nine 
AGN in the sample. Again, this is not a surprising result since 
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Fig. 9. Left Distribution of restframe equivalent widths of the sample of 170 candidates. Bin size is 10 A. The solid line indicates the 
exponential fit to the data with the fits with ±1<t plotted with dotted lines. The dashed line is the best fit to the sample of Gronwall 
et al. (2007). This sample has a more narrow distribution of equivalent widths. Right Distribution of restframe equivalent widths for 
total (non-GALEX-detected), red and blue sample (with Vj— + > or < 0). In this plot, only objects with detections in both the Vj 
and the i+ bands are given. Red sample is plotted with a red line and blue with a blue line. The width of the different distributions 
are given in the figure. 



redshift z 2 corresponds to the peak of the AGN number den- 
sity distribution (see e.g. Wolf et al. 2003 and references therein). 
The number density of AGN in this sample is in principle con- 
sistent with surveys at higher redshift, but is likely to be higher 
by a factor of two. 

- Change in SFR ratios 

The ratio of UV to Lya derived SFRs has a higher median value 
in this survey than results at redshift z ^ The spread is large, 
but the trend is clear. This result, and the smaller EWs discovered 
in general for this sample (see also the next point) at high redshift 
may indicate that LAEs in this survey are expected to be more 
affected by dust than LAEs at z ^ 3. 

- Narrower EW distribution 

Another intriguing and related result is that of the narrower 
equivalent width distribution. Firstly, as detailed in section 14.41 
the distribution is difficult to explain without invoking compli- 
cated arguments about the properties of the galaxies. A possible 
explanation would be that the star formation histories of Lya 
emitters are complex. Secondly, the difference in the distribu- 
tions between redshifts three and two is further evidence of a 
higher dust quantity in the lower redshift galaxies. 

In conclusion, by comparing observations of Lya emitters 
at redshift z = 2.25 with galaxies selected in the same manner 
at higher redshifts, several evolutionary signatures become evi- 
dent in the properties of these galaxies. At lower redshifts, there 
appear to be fewer objects, with redder colours and higher dust 
contents, smaller equivalent widths, and a higher fraction of ob- 
jects containing AGN. Future SED fitting of these galaxies will 
reveal more information into the properties such as dust, mass 
and age (Nilsson et al., in prep.). 
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